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1 Introduction 



One of the most interesting aspects of deep inelastic lepton-proton scattering is the 
contribution to the Qi spin-dependent structure function from photon-gluon fusion 
subprocesses '~i*{q)g{p) — > qq- Naively, one would expect zero contributions from 
light mass qq pairs to the first moment /q^ dx gi{x, Q'^) since the q and q have opposite 
heli cities. In fact, this is not the case if the quark mass is small compared to 
a scale set by the spacelike gluon virtuality p"^. This is the origin of the so-called 
anomalous correction —3^Ag to the Ellis- Jaffe sum rule for isospin-zero 

targets assuming three light flavors. Here Ag is the helicity carried by gluons in the 
hadron target, Ag{Q) = Jq dx[g'i{x,Q) — gi{x,Q)], at the factorization scale Q. In 
the language of the operator product expansion, the photon-gluon subprocess contri- 
butions to the first moment of gi{x,Q) correspond to the anomalous VVA triangle 
graph ip, 1^ contribution to the hadronic matrix element of the local axial current. 

For fixed gluon virtuality = —p^ the photon-gluon fusion process induces two 
distinct contributions to the first moment of gi in polarized deep inelastic scattering. 
Let nig denote the mass of the struck sea quark. When is much greater than both 

and the box graph contribution to the first moment of gi for a gluon target is 



2m2 1 Jl + AmyP^-1 

1 + -— , In 



^ 271 ^ ' P2 ^ 4ml/ P^ V \/l + 4m2/p2 + 

The first, mass- independent term (— ^) in Eq. ([l|) comes from the region of phase 
space where the struck quark carries large transverse momentum squared ~ 
relative to the photon-gluon direction. It measures a contact photon-gluon interaction 
and is associated with the axial anomaly |^. The second mass-dependent 
term comes from the region of phase space where the struck quark carries transverse 
momentum ~ P'^,ml. This mass dependent term vanishes in the limit P^ ^ 
and tends to +^ when ^ m^. The "soft" mass dependent term in Eq. (|1|) is 
associated with the quark parton distribution of the gluon Ag*^^^"°'^) ; it can safely be 
neglected for the light (up and down) quarks. 

On the other hand, the magnitude of the gluon virtuality is important for gauging 
the contribution of the massive sea quarks. If the sea quark mass is heavy compared 
to the gluon virtuality 4m^ ^ P^ = — p^, the photon-gluon fusion contribution 
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to dx gi{x,Q'^) vanishes to leading order in as{Q'^)- This result follows from a 
general theorem based on the Drell-Hearn-Gerasimov sum rule which states that 
the logarithmic integral over the photoabsorption cross section 



^ Aa^a^W^) = 0{a') ; (2) 



it vanishes at order for any 2 — > 2 Standard Model process ||T0|, |TT|. Here Act is 
the cross section difference for parallel versus anti-parallel incident helicities. In the 
present application, the gluon (for = 0) takes the role of the on-shell photon 7 
and the particle a can be taken as a real or virtual photon. As the photon virtuality 

becomes large, the DHG integral evolves to the first moment of the helicity- 
dependent structure function gi{x, Q"^). Thus the fusion ~^ Born contribution 
to dx gi{x,Q'^) vanishes for small gluon virtuality ^ 4m^, ^ Q"^. Notice 
that the Born photon-gluon fusion contribution to the Ellis- Jaffe moment is zero even 
for very light quarks as long as the gluon virtuality can be neglected. 

The above application of the DHG theorem holds for any photon virtuality = 



— Q , and is thus more general than leading twist [12|. In fact, the leading-order 



fusion contribution to the du/u moment of the difference of helicity-dependent photo- 
absorption cross sections vanishes even if < Am'^, as long as the gluon virtuality 
can be neglected. The result also holds for the weak as well as electromagnetic current 
probes P, IH. 

It is clearly important to ascertain the actual numerical contribution of heavy 
quarks ss, cc, bb, ti to the first moment of gi, i.e., what is the effective number of sea 
quark contributions to the Ellis- Jaffe moment? From the above discussion, the specific 
contribution of a given sea quark pair qq depends not only on Q^, but more critically 
on the ratio of scales p^/4m^. In a full QCD calculation of photon-gluon fusion 
contributions to the first moment of gi one needs to integrate over the distributions 
of extrinsic and intrinsic gluon virtualities in the target nucleon. For small gluon 
virtualities (P^ <^ m^) the "hard" anomaly contribution to the first moment of gj ^ 
cancels with the "soft" mass dependent contribution. For deeply-virtual gluons the 
mass-independent anomaly contribution dominates over the mass- dependent term 
which tends to zero. Therefore, we shall investigate the effect of retaining the finite 
quark masses and performing a more exact analysis, in which we integrate over P^. 
Our aim is to understand the role of heavy quarks (e.g. strange and charm) in the 
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photon-gluon fusion process. 

As we shall show in the next section, the exact form of the spectrum N{p^) of gluon 
virtuality in the target nucleon depends in detail on the physics of the nucleon wave- 
function. "Extrinsic" gluon contributions, which arise from gluon bremmstrahlung 
qv — > qvQ of a valence quark, have a relatively hard virtuality dNextip^) / dp^ ~ 
ois{p'^)/p'^, above a minimum virtuality p^in- '^^^ mean virtuality of the extrinsic 
gluons depends on the upper limit of integration, which in turn depends on the kine- 
matic phase space. On the other hand, intrinsic gluons, which are associated with 
the physics of the nucleon wavefunction (for example, gluons emitted by one va- 
lence quark and absorbed by another quark), have a relatively soft spectrum. We 
will characterize the shape of the intrinsic gluon virtuality by the convergent form 
dNintip^) / dp^ ~ dN extip^) / dp^ / [1 where is a typical hadronic mass 

scale. We shall use such model forms for the extrinsic and intrinsic gluon distribu- 
tions to predict specific contributions of the heavy sea quarks to the first moment of 
c/i(x,g2). 

In addition to the photon-gluon fusion contributions, additional contributions to 
the first moment of gi{x^Q'^) arise from intrinsic heavy sea quarks associated with 
higher Fock states in the target hadron. For example, meson-baryon fluctuations such 
as, p KA imply a negative intrinsic strange quark contribution to Jq dx gi{x,Q'^) 
[|14| . In the case of charm, the small probability 0(1%) of intrinsic charm present in 
the proton implies a small intrinsic charm contribution to Jq dx gi{x, Q"^). 

The charm contribution to the nucleon helicity-dependent structure functions and 
sum rules will be addressed by several new experiments. The COMPASS [jl5[ and 
HERMES [|I6| experiments will measure charm production |T^-[^ in polarized deep 
inelastic scattering. Experiments have also been proposed at SLAG The aim of 
these experiments is to learn about the gluon polarization in a nucleon through the 
photon-gluon fusion process. 



2 Polarized gluons and gi 

In order to analyze the sensitivity of the anomaly to the sea quark mass in the 
photon-gluon fusion subprocesses, let us start by expressing the contributing gluon 
distributions in terms of the corresponding bound-state wavef unctions. In general 
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the dependence of the parton distributions comes from the integral of the bound- 
state wavef unction over the virtuahty of the corresponding parton up to the scale Q^. 
Schematically, for the polarized gluon distribution, we have 

AG(x,Q2) = dP'[\^,y,,{P\xf-\^,y,^{P\x)\\ (3) 
which means that 

^ AG(x,P^) = \^,yAP\^)? - \^.ilAP\^)\' - (4) 



Qp2 



Here ^E'gt/p| and gi/p^ are the gluon wavefunctions for positive and negative helicities 
relative to the proton helicity as functions of the gluon virtuality = —p^ and the 
fraction x of the plus component of the target nucleon's momentum. 

In perturbative QCD the total photon-gluon fusion contribution to gi for a nucleon 
target is given by 

9Pi^,Q') =lT.4 9f'\x,Q'), (5) 

where g[^'^'^ is the contribution where the struck quark carries flavor q 
.io.)f^ n2, _ f^' .02 diAGix,P^))^ , 



gr'{x,Q') = J^^ dP'^^^^^^^(^A,{x,Q\P'). (6) 



min 



Here ® denotes the convolution over x and Aq denotes the contribution to the spin 
structure function gi of a "gluon target" with virtuality P^, where the struck quark 
carries flavor q. The infra-red cut-off -Pmin is the minimum gluon virtuality at which 
we can apply perturbative QCD — that is, where the current-quark and gluon degrees 
of freedom in perturbative QCD give way to dynamical chiral symmetry breaking and 
conflnement. The GRV and Bag model analyses of deep inelastic structure 
functions involve taking a QCD-inspired model-input for the leading twist parton 
distributions at some low scale /Xq, evolving the distributions to deep inelastic and 
comparing with data. The optimal GRV and Bag model flts to deep inelastic data 
are found with /Xq ~ 0.2 - 0.3 GeV^. Motivated by this phenomenological observation, 
we shall set P^.^ = 0.3GeV2. 

In the Born approximation Aq is calculated from the box graph contribution to 
photon-gluon fusion. We can deflne the "hard" part of Aq by imposing a cutoff on 
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the transverse momentum squared of the struck quark > 

/-I 4(m|^ 

Alw(x,Q^p^A^)- 



27T 1 



4^2p2 



{2x-l){l--^) (7) 



1 , /I 4x2 P2 4(m2+A2) 



1 In 



Here is the fermion mass, x is the Bjorken variable and W'^ = Q'^{^-^) — is 
the center of mass energy for the photon-gluon colhsion. The running couphng, as, 
in Eq. (7) is evaluated at the scale P^. Following Parisi and Petronzio we shall 
use a modified running as(P^) — see Eq. (15) below — which freezes in the infrared, 
to describe Aq when P^ becomes small. 

Keeping contributions where the struck quark carries transverse momentum squared 
> A^, the photon-gluon fusion contribution to the first moment of Qi^'^^ is obtained 
from (6): 

^,(g^ A^) = /J dP^ up\ A^) ^^(P^) . (8) 

min 

X,(P2, A^) = dx A(x, Q^ P^ A^) (9) 

JO 



Here 



and 



The cutoffs Xmax and Zmax in Eqs. (p|,p^ come from the kinematics. For the box 
graph term Aq, the cutoff Xmax = Q"^/ (Q^ + P^ + 4(mg + A^)) in Eq. (||) is obtained 
from the phase space factor \/l — — |/t^^^ in Eq. (H). The cutoff Zmax in Eq. (|iy) is 
derived from the explicit form of the polarized gluon distribution — see below. 

In the rest of this Section we discuss the contribution of qq pairs with small trans- 
verse momentum (when we relax the A^ cut-off), the size of higher-twist contributions 
to the first moment of gi^'^\ and the jet signature of the different contributions to 
the first moment. 

When we integrate over the full range of possible impact parameters we need to 
include small values of in Eqs. (P^TOl). This necessarily involves extrapolating 
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the calculation into the domain of non-perturbative QCD. Shore and Veneziano 
have considered the analogous process of the spin structure function of the polarized 
photon for a virtual photon target. They argue that the target photon virtuality 
where Jq dxgj grows from zero (at = 0) to — fA^c depends on the realization of 
chiral symmetry breaking in QCD. In perturbative QCD the individual quark flavor 
contributions to Jq dxgKx, Q^) grow rapidly from zero when ~ Am'^. In full QCD, 
spontaneous chiral symmetry breaking means that the scale of the transition virtuality 
is set by the constituent quark mass rather than by the current quark mass — that is, 



we expect Jq dxg1{x, Q"^) to grow rapidly from zero when P^ 



Motivated by this 



result, one might expect the gluon- virtuality where 1q grows rapidly to depend on 
any possible diquark structure of the gluon at low k"^. Spontaneous chiral symmetry 
breaking is a considerably more dramatic effect in the light quark masses than it is 
in the heavy quark masses. Thus we can expect that perturbative QCD will provide 
a reasonable model-independent estimate of the heavy-quark Aq when becomes 
small. 

When — > oo the expression for Aq simplifies to the leading twist (=2) contri- 
bution: 



Aq{x,Q\P\\^) 



as 
2^ 



(2a;-l)fln^ + ln^^-l' 



A2 



X 



(11) 



-{2x - l)ln 



X(l - X)P2 + (^2 ^ ^2 



+ (1 -xl 



2ml 



P^x{2x 



x(l-x)P2 + r7i2 + A2 



which has the first moment 

2m' 



27r 



1 + 



p2 



A2)/P2 



:ln 



+ 4(m2 + A2)/p2_ 
^l + 4(m2 + A2)/p2 + 1 



(12) 

For finite quark masses, the cutoff Xmax protects Aq from reaching the ln(l — x) 
singularity in Eq. (H). To quantify this effect, in Table |l| we list the values of Xq for 
different values of P^ and A^. The "Q^ = oo" values are obtained by keeping only 
the leading twist contribution, Eq. ([131) . For the "hard" cut-off A^ = 1 GeV^ the 
cut-off itself acts as a major source of higher-twist. When we relax the cut-off by 
setting A^ to zero we find a large ~ 63% higher twist suppression of at Q"^ = 10 
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Table 1: Heavy quark effects in Xq (in units of — "^^^ ' ) 



light strange charm 



0.5 1.0 10.0 0.77 0.74 0.16 

0.5 1.0 100.0 0.98 0.94 0.30 

0.5 1.0 oo 1.00 0.96 0.32 

0.5 0.0 10.0 0.99 0.61 0.013 

0.5 0.0 100.0 1.00 0.63 0.033 

0.5 0.0 oo 1.00 0.63 0.035 



GeV . Steffens and Thomas pj] have observed that Ic is, to good approximation, 
independent of for values of P^ between zero and 1 GeV^. The rise in Ic with 
increasing corresponds to removing a greater amount of the mass dependent term 
in Eq. (|l|) which cancels against the mass-independent (anomaly) term which is 
associated with ~ in the limit — oo. When we increase the cut-off on 
the transverse momentum squared of the struck quark we increase the infrared cut-off 
on the invariant mass of the qq pairs produced by the photon-gluon fusion. Quark 
mass dependent contributions to Xg go to zero when we increase the invariant mass 
Aiqq much greater than Arriq. Note that the light-quark X/ in Table ^ is significantly 
suppressed below unity with the A^ = 1 GeV^ cut-off at = 10 GeV^. If we 
decrease the cut-off this Ii grows to 0.87 (A^ = 0.5 GeV^), 0.91 (A^ = 0.3 GeV^) and 
0.96 (A^ = 0.1 GeV^) when we use our modified as ( — see Eq. (|T5|) below) together 
with the current light-quark mass through-out. 

Consider the large Q'^ limit (Q^ 3> 4m^). When we set A^ = in Eq. (|12D to obtain 
Eq. ([H), the "hard" anomaly contribution to Iq (the first moment of gj ^) cancels with 
the "soft" mass dependent contribution for small gluon virtualities (P^ ^ ^g)- For 
deeply virtual gluons the mass-independent anomaly contribution dominates over the 
mass- dependent term which tends to zero in the limit P^ ^ 4m^. It is interesting to 
observe that in a semi-inclusive experiment one can in principle identify events which 
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correspond specifically to the contributions to the first moment of gi{x). These are 
events with three jets recoiling, taking up the large momentum qt transferred by the 
lepton. The final state consists of the qq pair with perpendicular momentum qr — Pt 
plus the quark that emitted the gluon, with mass M and transverse momentum p^. 



These events have gluon virtuality —p > m^. As shown in |^T[, the events where 
only the qq pair recoils produce no contribution to the first moment of gi{x). This 
corresponds to events with small transverse momentum of the quark that emitted 
the gluon p"^ <^ m^, or with gluon virtuality —p'^ <^ m^. The vanishing of the first 
moment of gi{x) for heavy quark production implies that there must be a polarization 
asymmetry zero, which in principle can be measured experimentally pl[|. 



3 Extrinsic and intrinsic glue 

Having established the theoretical framework, we now investigate photon-gluon fusion 
using simple models for the exclusive and inclusive gluon distributions of the nucleon. 

The "extrinsic" glue consists of gluons which are radiatively generated from in- 
dividual valence quarks in the target whereas the "intrinsic" glue is associated with 
gluon exchange between valence quarks. For example, consider a gluon which is ex- 
changed between two valence quarks in the proton. In constituent quark models these 
"gluon exchange currents" contribute to the proton-Delta mass difference [^| . They 



also renormalize the valence contributions to the nucleon's axial charges, which are 
measured in /5— decays and in the first moment of gi{x, Q^). When cut, the exchanged 
gluon gives an intrinsic gluon. A gluon which contributes to the quark self energy 
when cut gives an extrinsic gluon. 

We now estimate the size of the extrinsic and intrinsic gluon contributions to the 
first moment of gi. We calculate the ratio Ues = ^hf^i of the heavy to light quark 
contributions to photon gluon fusion for both the extrinsic and intrinsic glue. 

The extrinsic and intrinsic gluon distributions are dominated by gluons with small 
virtuality. The virtuality distribution ^^^^ for the extrinsic glue contains a logarith- 
mic tale extending to the kinematic limits. Phenomenologically [^, the momentum 
distribution of the intrinsic glue is found to be weighted by the factor 

WiP') = NJ-^) (13) 
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relative to the extrinsic glue. The mass parameter can be estimated as 0.71 
GeV^, the mass scale of the dipole fit to the proton form factor [P3[; Ng is a model 
dependent normalization constant. 

We start with a simple model for the gluon distributions taking into account their 
correct distribution. Model dependent normalization uncertainties cancel in the 
ratios T^/Ti. We treat the nucleon target as a three-quark system where the target 
quarks are treated as "elementary" with constituent quark mass M equal to 300 MeV. 
The polarized extrinsic gluon distribution is given by 

P) - MC^- j + ^2^2)2 (1 - (14) 

where the {1 — z) factor is a Jacobian factor for the change of variables from to 
P2 in Eq. d). The QCD factor Cp = ^f ^ 0.6 H, H, H is the spin depolar- 
ization factor found in relativistic quark models, which parametrizes the transfer of 
the proton's angular momentum from intrinsic spin of the quarks to orbital angu- 
lar momentum through relativistic effects and quark-pion coupling. Since the gluon 
transverse momentum squared j9y = P^(l — z) — M'^z'^ is non-negative, we obtain the 
W cutoff in Eq. (|10D: z^^,{P^) = (-1 + Vl + 4C)/2C, where C = M^P^ 

Our simple model, Eq. (|14D , for the gluon distributions exhibits the a; 
behavior predicted by color coherence By construction, it also exhibits the large 
X behavior associated with an elementary quark target. In a more sophisticated model 
one should also include gluon exchange between the valence quarks in addition to the 
gluon involved in the qq process. However, this is beyond the scope of the 

present paper. 

We work in the analytically extended |]38| ay scheme This means that we 
use the running coupling 

Att 

^siP') = . ,p.^w. . (15) 



in Eqs. (7,14). Here = 0.2GeV^, Ay = 0.16GeV and the number of flavors 
which contribute to (3q is taken as a continuous variable which depends on P^ 



/3o = (11 - I E N,) where 
i=i 



N,^[l + j^ \ (p, = P^/ml) . (16) 
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Figure 1: The effective number of flavors ngfr for heavy sea quarks ss, cc, and bb 
contributing to the first moment of gi{x, Q"^), arising from 7*-(extrinsic gluon) fusion, 
as a function of momentum transfers < lOOGeV^. In Fig. |l]a, the cutoff on quark 
transverse momentum > is set equal to zero. In Figs. |l]b and |l|c, =1 GeV^ 
and A^ = 10 GeV^, respectively. 
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Table 2: The effective number of heavy-flavors Ues = ^h/^i- 
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In Figs. |T]a-|T]c we show the effective number of flavors, ries = ^hf^h for the heavy 
flavor ss, cc and bb production contributions to the first moment of gi{x^Q'^) from 
7*-(extrinsic gluon) fusion up to = lOOGeV^. Figure ^ is obtained by setting 
equal to zero. (In this calculation we have used "modified" as, Eq. (^^, together 
with the current light-quark mass through-out.) In Figs. and |l]c we set the cut-off 
A^ equal to 1 GeV^ and 10 GeV^ respectively. We repeat these calculations for 7*- 
(intrinsic gluon) fusion in Figs.0a-^. The results in Figs. |1] and |^ are summarized 
in Table |^. 

The effective number of flavors Ucs = F/i/F; increases for the heavy quarks when 
we increase the cutoff A^ on the transverse momentum squared of the struck quark. 
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Figure 2: The effective number of flavors Ues for lieavy sea quarks ss, cc, and bb 
contributing to the first moment of gi{x, Q^), arising from 7*-(intrinsic gluon) fusion. 
In Fig. I^a, the cutoff on quark tranverse momentum k'^ > is set equal to zero. In 
Figs. |b and |c, =1 GeV^ and = 10 GeV^, respectively. 
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Figure 3: The integrand in Eq. as a function of for the hght and charm quark 
contributions at = 100 GeV^ for extrinsic glue. Figure ^3 shows the corresponding 
integrand for intrinsic glue. 

This corresponds to removing a greater amount of the mass dependent term in Eq. 
([1|) which cancels against the mass-independent (anomaly) term from ~ in the 
limit Q"^ —>■ oo. By increasing we are increasing the cut on the invariant mass of the 
qq pairs produced by the photon-gluon fusion. Quark-mass dependent terms become 
less important when the invariant mass Aiqq becomes much greater than 4m^. 

In Fig. ^ we show the integrand in Eq. (|^) as a function of for the light 
and charm quark contributions at = 100 GeV^ with the extrinsic glue. Figure 
^3 shows the corresponding integrand for the intrinsic glue. Figures ^ and both 
involve = 0. Note that the heavy and light quark curves come closer together with 
increasing P^. This result corresponds to the fact that the mass-dependent term in 
Eq. (|I]) tends to zero in the limit P^ ^ Arri^. 
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4 Phenomenology and discussion 



Gluon polarization offers a possible explanation for the small value of g^^^ (the three- 



flavor, singlet axial charge) extracted from polarized deep inelastic scattering ^ 



~ 0.2 - 0.35. (17) 



Relativistic binding |Q and constituent-quark pion coupling [0, Q models predict 
(yf^-* ~ 0.6 — a factor of two larger than the measured g^^\ In these semi-classical 
models g^^^ is interpreted as the fraction of the nucleon's helicity which is carried 
by its quark constituents. In QCD the axial anomaly ||^, ^ induces various gluonic 
contributions to g^^\ One finds |jl|, 0, H, ^| 



+ c . 

partons 



Here ^Aq and Ag are the amount of spin carried by quark and gluon partons in the 
polarized proton. The — S^A^fpartons term is associated with the mass- independent, 
local 7*5' interaction in Eq. (|1|) assuming three light flavors. The soft mass-dependent 
contributions to photon-gluon fusion are included in Agpartons- The last term, C, 
is associated with non-trivial gluon topology and a possible 6{x) term in gi. 
It is missed by polarized deep inelastic scattering experiments which measure the 
combination {g^^^ — C). 

How large are the photon-gluon fusion sea-quark contributions Tg to g^^^ if we 
allow for finite sea-quark masses and a spectrum of gluon virtuality? 

Since extrinsic glue is radiatively generated from single quark lines in the target, 
we believe that our model should provide a good order-of-magnitude estimate for the 
normalization of the extrinsic Tq. We find * = -0.0024 and (r„ + Td + T,)^^* = 
-0.033 at = 100 GeV^ for the extrinsic glue. 

The magnitude of the intrinsic gluon contribution to g^^^ depends on the normal- 
ization of the gluon polarization. Taking the estimate A(yfintrinsic = +0.5 at 1 GeV^ 
13, we obtain rj°* = -0.0020 and (r„ + + r,)'°* = -0.07 when = loo GeV^. 

It is interesting to compare these estimates of with the results of heavy-quark 



effective theory |21] and heavy-quark operator product expansion calculations. 
These calculations express the total heavy-quark contribution to the first moment of 
gi in terms of the three-light- flavor singlet axial-charge g^^\ Using the most recent 
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value, Eq. (|T^, of g)^ , Manohar's effective theory calculation of the heavy-charm- 
quark axial-charge becomes 

^(charm)^g2 ^> ^ -0.0055 ± 0.0018 + 0(1/771^) . (19) 

Manohar gives an estimate ~ 0.003 (magnitude) for the 0{l/mc) corrections. 

The sum of our extrinsic and intrinsic charm contributions (r°^* + r™* = —0.0044) 
is in good agreement with Eq. (|1^). However, with the same gluonic input, photon- 
gluon fusion can account for only about one-third of the difference between the value of 
g^^^ extracted from polarized deep inelastic scattering and the quark model prediction. 
Next-to- leading order QCD fits to the present world data for gi are consistent with a 
value of Ag between zero and +2 at 1 GeV^ [^. The value Ag = 2 would increase 
our estimate of the intrinsic gluon contribution by a factor of 4 and would bring 
theory into agreement with the empirical determinations of g^^\ We look forward 
to a more precise measurement of Ag from forthcoming experiments on open charm 
production. 

Finally, it is interesting to note that since the contributions due to heavy sea 
quarks come from highly virtual gluons, one expects minimal nuclear shadowing for 
their contribution to the first moment of gi . 
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